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The  first  two  dimensional  mathematical  model  of  a  complete  direct  propane  fuel  cell  (DPFC)  is  described. 
The  governing  equations  were  solved  using  FreeFem  software  that  uses  finite  element  methods.  Robin 
boundary  conditions  were  used  to  couple  the  anode,  membrane,  and  cathode  sub-domains  successfully. 
The  model  showed  that  a  polytetrafluoroethylene  membrane  having  its  pores  filled  with  zirconium  phos¬ 
phate  (ZrP-PTFE),  in  a  DPFC  at  150°C  performed  much  the  same  as  other  electrolytes;  Nation,  aqueous 
H3  PO4,  and  H2  SO4  doped  polybenzimidazole,  when  they  were  used  in  DPFCs.  One  advantage  of  a  ZrP-PTFE 
at  150°C  is  that  it  operates  without  liquid  phase  water.  As  a  result  corrosion  will  be  much  less  severe 
and  it  may  be  possible  for  non-precious  metal  catalysts  to  be  used.  Computational  results  showed  that 
the  thickness  of  the  catalyst  layer  could  be  increased  sufficiently  so  that  the  pressure  drop  between  the 
reactant  and  product  channels  of  the  interdigitated  flow  fields  is  small.  By  increasing  the  width  of  the 
land  and  therefore  the  reactant’s  contact  time  with  the  catalyst  it  was  possible  to  approach  1 00%  propane 
conversion.  Therefore  fuel  cell  operation  with  a  minimum  concentration  of  propane  in  the  product  stream 
should  be  possible.  Finally  computations  of  the  electrical  potential  in  the  ZrP  phase,  the  electron  flux  in 
the  Pt/C  phase,  and  the  overpotential  in  both  the  anode  and  cathode  catalyst  layers  showed  that  serious 
errors  in  the  model  occurred  because  proton  diffusion,  caused  by  the  proton  concentration  gradient,  was 
neglected  in  the  equation  for  the  conservation  of  protons. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Several  advantages  are  obtained  when  hydrocarbons  are  the  fuel 
used  directly  at  the  anode  of  a  fuel  cell.  The  capital  cost  for  process¬ 
ing  equipment  needed  to  produce  hydrogen  or  methanol  can  be 
eliminated.  The  infrastructure  to  deliver  hydrocarbons  is  already  in 
place  (natural  gas  pipelines  in  urban  areas  and  commercial  trucking 
systems  for  delivery  of  LPG  and  diesel  fuel  in  rural  areas).  A  substan¬ 
tial  additional  capital  cost  for  delivery  systems  would  be  required  if 
hydrogen  and/or  methanol  were  to  become  equally  accessible.  Fur¬ 
thermore,  the  conversion  energy  inefficiencies  (e.g.  endothermic 
heat  required  for  steam  reforming)  are  eliminated.  Storage  of  liquid 
hydrocarbons  is  also  much  more  energy  efficient  and  convenient 
than  that  of  hydrogen. 

Although  hydrocarbons  have  many  advantages  they  also  have 
one  major  disadvantage  compared  to  hydrogen  or  methanol. 
Hydrocarbon  reaction  rates  (current  densities)  in  fuel  cells  are  more 
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than  one  order  of  magnitude  slower  than  those  of  hydrogen  or 
methanol. 

Extensive  fuel  cell  research,  using  hydrocarbons  as  the  fuel 
was  performed  in  the  1959-1968  period.  That  work  has  been 
summarized  in  three  reviews,  Liebhafsky  and  Cairns  [1],  Bockris 
and  Srinivasan  [2],  and  Cairns  [3].  Different  hydrocarbons  such 
as  normal  paraffins  (methane  to  n-hexadecane),  olefins  (ethylene, 
propylene  and  1 -butene),  gasoline  and  diesel  have  been  investi¬ 
gated  as  the  fuel  cell  feed.  An  overview  [1]  of  work  during  the 
1959-1968  period  shows  that:  (1)  normal  alkanes  have  been  the 
most  frequently  used  fuels,  (2)  most  electrolytes  were  aqueous 
acids,  (3)  platinum  has  been  the  favourite  electrocatalyst,  and  (4) 
cell  potentials  have  been  low,  which  is  a  serious  drawback. 

Work  in  our  laboratory  is  directed  toward  improving  the  viabil¬ 
ity  of  direct  hydrocarbon  fuel  cells  by  combining  several  strategies. 

The  first  is  the  selection  of  a  membrane  that  is  capable  of 
operation  at  a  temperature  of  at  least  150°C  to  ensure  that  no 
water  is  present  in  the  liquid  phase.  The  membrane  is  composed 
of  two  components,  a  porous  polytetrafluoroethylene  (PTFE)  that 
contains  zirconium  phosphate  (ZrP)  in  its  pores.  It  was  chosen 
for  this  research  because  it  has  demonstrated  reasonable  perfor¬ 
mance  when  operated  at  a  120  °C  in  a  hydrogen  fuel  cell  [4].  In 
the  absence  of  liquid  water,  corrosion  problems  are  expected  to 
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Nomenclature 

Apt  platinum  surface  area  per  catalyst  volume 

(mPtmc“a3talyst) 

Azat  specific  surface  area  of  catalyst  support  (Vulcan  car- 
bon)  in  the  anode  and  cathode  (m2ata|yst  g"1^) 
cG  total-concentration  of  the  gas  phase  (mol  m-3 ) 

Q, ely  concentration  of  charged-species  i  in  electrolyte 

(molm-3) 

Dj  diffusion  coefficient  of  species  i  (cm2  s-1 ) 

Dp  effective  particle  diameter  (fxm) 

F  Faraday’s  constant,  96485  (C  mol^^) 

/pt/CAT  fraction  of  catalyst  support  surface  area  covered  by 
Platinum  (m2tm-2talyst) 

AG*  activation  energy  for  the  exchange  current  density 
(kjmol-1) 

j  volumetric  current  density;  rate  of  production  of 

proton  in  electrode  (Am-3talyst) 

j  current  density,  vector  quantity  (A  m-2ctrode  face  area ) 

jo  exchange  current  density  at  operating  conditions 

(Amp2) 

Jo-ref  reference  exchange  current  density  at  the  reference 
conditions  (Am“2) 

J  total  current  density  ( mA  cm-2 ) 

Ji  diffusion  flux  of  species  i  (mol  m-2  s-1 ) 

Lw  land  width  in  the  flow  field  (mm) 

MW/  molecular  weight  of  species  i  (g  mol-1 ) 

n  normal  direction 

N/  ely  flux  of  species  i  in  electrolyte,  vector  quantity 
(mol  cm-2  s-1) 

Pi  partial  pressure  of  species  i  (Pa) 

P  total  pressure  (Pa) 

R  universal  gas  constant  (J  mol-1  K-1 ) 

Si  consumption  rate  of  the  reactant  species  i 

(molm-3  s-1) 

T  temperature  (K) 

Th  thickness  of  catalyst  layers  and  membrane  (p,m) 

t  time  (s) 

u  velocity  of  gas  mixture,  vector  quantity  (m  s-1 ) 

uEly  velocity  of  electrolyte  (ms-1) 

U[  mobility  of  species  i  (cm2  molj-1  s-1 ) 

un  velocity  of  gas  mixture  in  normal  direction  (m  s-1 ) 
x  Cartesian  coordinate 

y  Cartesian  coordinate 

yi  mole  fraction  of  species  i  in  the  gas  phase 

z  moles  of  electron  in  anode  and  cathode  reactions 

(moleiectrons  Uiolpropane) 

Zj  charge  number  of  species  i  (molcharge  mol“pecies) 

Greek  letters 

aA  and  ac  anodic  and  cathodic  charge  transfer  coefficients 
,  ot2  and  a3  relaxation  parameters  in  boundary  conditions 
r]  overpotential  (V) 

s  volume  fraction 

p  mass  density  (kg  m-3 ) 

Azat  apparent  bulk  density  of  catalyst  support 

(Scatalyst  mCataIyst) 

Pi  mass  density  of  species  i  (kg  m-3 ) 

jiG  dynamic  viscosity  of  gas  mixture  (kg  m-1  s-1 ) 

Vj  stoichiometric  coefficient  of  species  i 

a  ionic  conductivity  (S  m-1 ) 

O  electrical  potential  (V) 


^Pt-EQ 

equilibrium  potential  of  conductor  phase  (V) 

^ELY-EQ 

equilibrium  potential  of  electrolyte  phase  (V) 

A^cell 

cell  potential  difference  (V) 

EM 

sum  of  the  atomic  volumes  for  each  of  the  i  species 
(m3  kg  atom-1) 

Subscripts  and  superscripts 

A 

anode 

av 

phase-averaged  quantity 

C 

cathode 

c3 

propane 

C3Ox 

propane  oxidation  reaction  on  Pt  catalyst 

C02 

carbon  dioxide 

ELY 

electrolyte  phase  in  the  anode  or  cathode  catalyst 
layers 

EQ 

equilibrium  state 

G 

gas  mixture 

h2o 

water 

i 

species  in  gas  phase;  propane,  water,  C02  and  02 

ML 

membrane  layer 

n 

normal  direction 

o2 

oxygen 

02Rd 

oxygen  reduction  reaction  on  platinum  catalyst 

Pt 

platinum  catalyst 

ref 

reference  conditions 

decrease  dramatically  and  it  may  be  possible  to  replace  precious 
metal  catalysts  with  less  expensive  ones.  The  model  described 
here  uses  Pt/C  catalysts  for  both  anode  and  cathode  because 
reaction  rate  data  were  available  in  the  literature.  If  the  ZrP-PTFE 
membrane  investigated  here  is  shown  to  be  suitable,  then  a  future 
search  for  non-precious  metal  catalysts  maybe  appropriate. 

The  use  of  interdigitated  flow  fields  has  at  least  two  desirable 
features  [5,6].  First,  the  partial  pressure  of  the  reactants  all  along 
the  feedstock  channels  can  be  maintained  constant  at  their  val¬ 
ues  in  the  feedstock.  In  the  more  common  serpentine  flow-fields, 
reactants  are  diluted  by  the  products  that  enter  the  same  channel. 
Therefore,  the  partial  pressures  of  the  reactants  decrease  as  reac¬ 
tion  proceeds  and  more  products  are  produced.  Second,  the  width 
of  the  catalyst  layer  can  be  increased  thereby  increasing  the  reactor 
length  and  the  residence  time  in  contact  with  the  catalyst.  In  prin¬ 
ciple,  the  conversion  could  approach  100%,  thereby  minimizing  the 
concentration  of  unreacted  hydrocarbons  in  the  product  stream. 

This  work  describes  the  first  two-dimensional  model  of  a  com¬ 
plete  direct  propane  fuel  cell  (DPFC).  Both  anode  and  catalyst  layers 
contain  two  solid  phases,  a  Pt/C  phase  that  conducts  electrons 
and  a  ZrP  phase  that  conducts  protons.  The  ZrP-PTFE  membrane 
conducts  protons  through  the  ZrP  that  fills  the  pores.  There  are 
no  other  components  in  the  membrane  layer.  The  following  reac¬ 
tions  proceed  in  a  DPFC  assuming  complete  propane  conversion. 
The  electrical  potential  values  shown  are  the  standard  reversible 
electrochemical  potentials  at  25  °C: 

Anode  reaction  :  C3FI8  +  6H20  =  3C02+20H++20e-  </>AN  =  0.136V 
Cathode  reaction  :  502  +  20H+  +  20e-  =  1OH2O0Ca  =  1.229V 

Overall  reaction  :  C3FI8  +  502  =  3C02+4H20  A0C ell=  1.093  V 

The  generation  of  electrical  energy  in  rural  areas  is  our  target 
application  for  DPFCs.  The  cost  of  delivering  electrical  energy  to 
rural  areas  is  substantially  greater  than  to  urban  areas,  because 
longer  transmission  lines  are  required  to  serve  a  comparatively 


3188 


H.  Khakdaman  et  al.  /  Journal  of  Power  Sources  196(2011)  3186-3194 


Anode  bi-polar  plate 
Anode  catalyst  layer 
Membrane  layer 
Cathode  catalyst  layer 


Cathode  bi-polar  plate 


Fig.  1.  A  direct  propane  fuel  cell  with  interdigitated  flow  field. 


small  number  of  customers.  Therefore  more  costly  fuel  cells  can 
be  justified  for  rural  areas  than  for  urban  areas.  In  addition,  the 
infrastructure  to  deliver  liquefied  petroleum  gas  (LPG)  or  propane 
to  rural  areas  already  exists,  because  LPG  is  one  of  the  fuels  used 
for  heating  in  rural  areas. 

2.  Model  development 

A  schematic  of  a  typical  DPFC  has  been  illustrated  in  Fig.  1.  The 
cell  is  composed  of  two  bi-polar  plates,  two  catalyst  layers  and  a 
membrane  layer.  The  cathode  bi-polar  plate  in  Fig.  1  shows  two 
sets  of  channels,  one  for  reactants  and  one  for  products  that  are 
connected  to  each  other  through  the  catalyst  layer.  In  addition 
the  catalyst  layers  can  be  thicker  when  using  interdigitated  flow 
fields  than  when  using  serpentine  flow  fields.  Thicker  catalyst  lay¬ 
ers  make  it  possible  to  have  a  better  dispersion  of  catalyst  (more 
reaction  sites). 

A  cross  sectional  view  of  the  cell  is  shown  in  Fig.  2.  Because  of 
the  symmetry  of  the  interdigitated  flow  field,  the  modeling  domain 
is  defined  to  be  from  the  middle  of  a  feed  channel  to  the  middle  of 
its  adjacent  product  channel.  It  contains  the  membrane  electrode 
assembly  (MEA)  including  both  the  anode  and  cathode  catalyst 
layers  and  the  membrane  layer. 

2. 1.  Governing  equations 

Anode  and  cathode  catalyst  layers  are  composed  of  three 
phases:  reactants  and  products  in  the  gas  phase,  solid  ZrP  elec- 


Propane  C02,  Water 

&  Water  &  Propane 


Air 


Water 
&  Air 


Fig.  2.  Modeling  domain  containing  anode  and  cathode  catalyst  layers  and  elec¬ 
trolyte  layer. 


Table  1 

List  of  equations  for  the  anode. 


Equations 


Conservation  of  mass: 

*S^=-V-(e  GPGU)-^ 
Conservation  of  momentum: 


6MWh2oJA  ,  3MWCo2  jA 

7F  +  7C 


-  VP  =  150 


mgU-sg  ) 

Dlel 


Conservation  of  non-charged  species: 

Gas  phase: 

—i fC3>  =  -  V  ■  (eGccuyc3  )+V-(Jc3)-j£ 

8(ecYH20)  =  “v  '  («GCGuyH2o)  +  V  •  (Jh2o)  - 

Conservation  of  charge: 

For  ionic  current :  V  ■  (sely^ely  V  0Ely,a)  =  -j a 
Butler-Volmer  equation: 

Ja  =  Jo.A-Apt  [exp  (  )  -  exp  (  )  ] 

where, 


,■  A  _  ,-c3ox  (  PC 3  A 
JO  A  J  o_ref,  \  PC3,ref  / 


exp 


rj A-  A0 A  -  A0EQ.A  -  (^Pt,A  -  ^ELY,a)  -  (^Pt-EQA  -  ^ELY-EQ.) 


Eq.  no. 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 


trolyte,  and  solid  catalyst  (Pt/C).  The  membrane  layer  also  contains 
solid  ZrP  as  well  as  PTFE.  Because  different  variables  in  different 
phases  are  of  interest,  a  multi-fluid  volume  averaging  method  has 
been  used  [8].  In  this  method,  spatial  averaging  is  performed  for 
each  single  phase  within  a  multiphase  control  volume.  Also,  the 
extensive  phase-averaged  quantities  are  related  to  the  intensive 
phase-averaged  amounts  and  phase  fractions.  For  example,  the 
phase-averaged  density  for  the  gas  phase  in  a  control  volume  can 
be  calculated  by  multiplying  the  gas  phase  density  by  its  volume 
fraction  in  the  control  volume,  as  written  in  Eq.  (1 ): 


PG,av  —  x  Pc  (1) 

Therefore,  phase-averaged  quantities  were  used  in  the  govern¬ 
ing  equations  for  the  solid  and  gas  phases.  Conservation  of  mass, 
momentum,  and  species  for  the  gas  phases  in  the  anode  and  cath¬ 
ode  are  shown  in  Tables  1  and  2,  respectively.  These  tables  also 
include  equations  for  the  conservation  of  protons  in  the  ZrP  solid 
phase  in  each  of  the  layers.  The  conservation  of  protons  is  the  only 
equation  that  has  to  be  solved  for  the  membrane  layer.  Because 
of  the  relatively  large  electron  conductivity  of  the  Pt/C  solid  cata¬ 
lyst  (the  electron  conducting  phase),  conservation  of  electrons  in 
the  Pt/C  phase  was  not  considered.  Instead  it  is  assumed  that  the 
electrical  potential  in  the  Pt/C  phase  was  constant. 


Table  2 

List  of  equations  for  the  cathode. 


Equations 


Conservation  of  mass: 

=-V.(£gPgu)- 
Conservation  of  momentum: 


5MWo2jc  |  10MWh2o/c 

7F  T  7F 


-  VP  =  150 


McO-ec) 


Conservation  of  species: 

Gas  phase: 

3(eGcGY02 ) 

—3 T2’  =  -V  ■  (sGcGiiyN2 )  +  V  •  Un2  ) 

Conservation  of  charge: 

For  ionic  current :  V  ■  (sely^ely  V  0Ely,c)  =  -jc 
Butler-Volmer  equation: 

jc=jo,c*P«[exp(^)-exp(^)] 

where, 


Jo,c=J 


;02  Rd  f  P°2  A 

'0-ref,  l  Po2,ref  J 


exp 


AG* 

02Rd 


TJc  =  A0C~  A0EQ.C  =  (^Pt,C  -  ^ELY.c)  -  (^Pt-EQC  ~  ^ELY-EQ.) 


Eq.  no. 
(10) 
(11) 

(12) 

(13) 

(14) 

(15) 

(16) 
(17) 
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Table  3 

Parameters  for  Eqs.  (8)  and  (16)  [10]. 


Parameter 

Anode 

Cathode 

Tref(K) 

423 

298 

Pref(kPa) 

101.3 

101.3 

Jo-ref  (A  Cm”2  Pt) 

1.0  xlO“8 

3.8  x  10-13 

AG*  (kjmol-1) 

90 

92 

The  first  term  of  the  Ergun  equation  has  been  used  to  describe 
momentum  conservation  because  the  catalyst  layers  are  essen¬ 
tially  packed  beds.  Conservation  of  non-charged  species  in  the  gas 
phase  accounts  for  diffusion,  convection  and  reaction  phenomena. 
The  consumption  rate  of  the  reactant  species,  S,-,  is  related  to  the 
volumetric  current  density,  j,  through  Faraday’s  law: 


St  = 


zF 


with  j  given  by  the  Butler-Volmer  equation: 


j  =  jo^pt 


(18) 


(19) 


where  jo,  the  exchange  current  density,  is  a  function  of  the  refer¬ 
ence  exchange  current  density,  jo-ref,  and  operating  temperature 
and  pressure.  The  reference  exchange  current  density  for  propane 
oxidation  on  platinum  has  been  reported  in  the  literature  [9].  In 
order  to  calculate  the  exchange  current  density  at  different  temper¬ 
atures  the  Arrhenius-type  expression  suggested  by  Psofogiannakis 
et  al.  [10]  was  used: 


Jo  =  Jo-ref  exp 


(20) 


The  parameters  required  to  calculate  j0  for  propane  oxidation 
at  the  anode  and  oxygen  reduction  at  the  cathode  are  reported  in 
Table  3.  The  final  expressions  for  the  exchange  current  densities  at 
the  anode  and  cathode  are  shown  in  Eqs.  (8)  and  (16). 

Three  mechanisms,  diffusion,  migration,  and  convection,  are 
involved  in  the  transport  of  protons  through  the  electrolyte  phase 
as  shown  by  Nernst-Planck  equation  [11]: 


Nu  ely  =  -ZjUjfCj  ELY  V  <2>ely  -  Dt  V q  ELY  +  q,  elyuely  (21) 

The  first  term  in  Eq.  (21)  describes  the  migration  of  protons  for 
which  the  driving  force  is  the  electrical  potential  gradient.  The 
second  term  describes  the  diffusive  flux  of  protons  for  which  the 
driving  force  is  the  concentration  gradient  (related  to  the  chemical 
potential  gradient).  The  final  term  describes  the  convective  flux  of 
protons  that  is  caused  by  the  bulk  motion  of  the  solvent.  For  the  sta¬ 
tionary  membrane  of  the  present  system,  uEly  =  0.  Some  of  the  more 
rigorous  models  have  included  both  migration  and  diffusion  terms 
[12-14].  However,  the  majority  of  PEM  fuel  cell  models  neglect  the 
concentration  gradient  and  therefore  the  diffusion  term  [  1 5-1 7].  In 
other  words  the  majority  of  PEM  fuel  cell  models  combine  migra¬ 
tion  and  diffusion  into  a  single  term  for  which  the  driving  force  is 
the  electrical  potential  gradient.  This  work  also  ignored  the  con¬ 
centration  gradient  and  described  proton  flux  as  a  single  term  for 
which  the  electrical  potential  gradient  was  the  driving  force.  That 
choice  had  some  undesirable  consequences  that  will  be  discussed 
in  a  subsequent  section  of  this  communication. 

Movement  of  charged  species  in  the  ZrP  electrolyte  phase  pro¬ 
duces  current  that  can  be  quantified  by  Eq.  (22)  [18]: 


CO 

3 

3 

CD 

3 


Anode  in  Wall  Anode  out 


Cathode  in 


Wall 


;  Cathode  out 


Fig.  3.  Boundaries  in  the  modeling  domain. 


where 

0 ELY  =  F2  zf  Ui q  ELY  ( 24 ) 


is  the  conductivity  of  the  ZrP  electrolyte. 

Electro-neutrality  of  the  entire  domain  can  be  expressed  by  Eq. 
(25): 


£v-*=° 


(25) 


As  the  charged  species  are  electrons  and  protons,  Eq.  (25)  at  any 
location  in  the  domain  can  be  written  as: 


-v  JH+  =  V  je_  =j  (26) 

where  j  is  rate  of  reaction  (volumetric  current  density)  given  by  the 
Butler-Volmer  equation.  Therefore,  an  equation  for  the  conserva¬ 
tion  of  protons  in  the  absence  of  electrochemical  reactions,  which 
is  the  situation  in  the  membrane  layer,  can  be  written  as: 


V  •  (<xelyV <Pely,ml)  =  0  (27) 

The  combination  of  Eqs.  (23)  and  (26)  is  shown  for  the  anode  as 
Eq.  (6)  and  for  the  cathode  as  Eq.  (14). 

2.2.  Boundary  conditions 

Fig.  3  shows  the  four  types  of  boundary  conditions  that  were 
used  to  model  the  domain;  inlets,  outlets,  walls  of  the  lands,  and 
symmetry  boundaries.  Table  4  shows  the  values  for  the  boundary 
conditions.  The  Ergun  equation  was  solved  using  the  boundary  con¬ 
dition  values  in  Table  4  for  the  velocities  at  the  inlets  (calculated 
from  inlet  flow  rates)  and  the  atmospheric  pressure  boundary  con¬ 
dition  at  the  outlets.  At  the  interface  between  the  catalyst  layers 
and  lands  (referred  to  as  the  wall),  the  flux  of  gaseous  species  is 
zero.  That  is  also  true  for  the  proton  flux  because  protons  are  not 
transferred  from  the  catalyst  layers  to  the  lands.  Instead  a  potential 
difference  is  set  between  lands  for  electrons.  The  zero  flux  condi¬ 
tion  has  been  applied  at  the  symmetry  boundaries,  assuming  the 
fuel  cell  repeats  by  reflection  through  these  boundaries  at  the  left 
and  right  of  the  computational  domain.  This  reflection  is  a  natural 
assumption  for  interdigitated  channels. 

The  feed  composition  is  known  at  the  inlet  of  the  catalyst  layers. 
It  is  assumed  that  no  change  in  the  composition  of  gaseous  species 
will  occur  after  the  gas  mixture  leaves  the  catalyst  bed.  Therefore, 
the  gradient  of  the  composition  is  zero  in  the  direction  normal  to  the 
catalyst  surface  at  the  outlet  boundaries.  Also,  the  electrical  current 
is  zero  in  the  normal  direction  of  the  inlet  and  outlet  boundaries, 
both  for  protons  and  electrons. 


j  =  FX’jZjNjely  (22) 

Substitution  of  Eq.  (21 )  in  Eq.  (22)  for  a  stationary  membrane  in 
the  absence  of  a  concentration  gradient  results  in: 

(23) 


2.3.  Model  input  parameters 

The  parameters  used  for  the  simulations  are  shown  in  Table  5. 
To  calculate  platinum  surface  area  per  catalyst  volume,  APt,  that  is 
used  in  the  Butler-Volmer  equation  it  is  assumed  that  a  fraction, 


j  =  -CTELyV^eLY 
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Table  4 

List  of  boundary  conditions  for  the  modeling  domain. 

Anode  inlet: 
uy  =  0.001  -0.005  m  s_1 
yC3  =0.10 
Anode  outlet: 

P=  101.3  kPa 
Cathode  inlet: 
uy  =  0.001  -0.005  m  s_1 
yo2  =  1.00 
Cathode  outlet: 

P=  101.3  kPa 

Anode/land  interface  (wall): 
uy  =  0 

Cathode/land  interface  (wall): 
uy  =  0 

Zero  flux  boundaries: 
ux  =  0 

Anode/membrane  interface: 

SELY.A^ELY.A  (  d0fy’A  )  =tfML  (^L) 

Un  =  0 

Cathode/membrane  interface: 

CTML  (  ^  )  =  £-ELY,C O’ELY, C  (  ) 

Un=0 


9^ELY  _  n 
dy  “U 

yw2  0  =  0.90 

9<%ly  _  n 

9y  “u 

5^1-0 

Uy  -u 

9<%ly  _  n 
dy  “U 

yn2o  —  0.00 

9<%ly  _  n 

9y  -u 

0 

II 

ol  >> 

TO1  1 

9<%ly  _  n 

gy  -u 

0 

II 

WUN  =  o 
dy 

0 

II 

9<%LY  _  n 
dx  U 

1=° 

§=o 

9Yh2o 


=  o 


0pt  =  0.1 5-0.50  V 


0Pt  =  O.9-1.2V 


dy± 

dx 


0 


^ELY.A  “  ^ML 


^ELY.C  =  ^ML 


/pt/CAT>  of  catalyst  support  surface  area,  ACat,  is  covered  by  platinum. 
Therefore  APt  can  be  calculated  using  Eq.  (28): 

Apt  —  /pt/CAT  x  ^CAT  x  PCAT  (28) 

2.4.  Numerical  procedure 

The  software  used  to  solve  the  two-dimensional  partial  dif¬ 
ferential  equations  is  FreeFEM++,  which  was  developed  by  Hecht 
et  al.  [19].  It  is  a  open-source  software  based  on  the  Finite 
Element  Method  which  is  capable  of  handling  multi-variables, 
multi-equations,  two  and  three-dimensional  systems  and  steady 
state  or  time  dependent  problems.  In  addition,  the  results  calcu¬ 
lated  using  FreeFEM++,  can  be  easily  exported  to  ParaView  software 
[20]  for  post-processing.  ParaView  is  open-source  software  used  for 
visualization  that  is  quite  powerful. 


Table  5 

Operational,  electrochemical  and  design  parameters  for  simulations. 


Property 

Value 

Temperature,  T 

423  K 

Pressure,  P 

101.3  kPa 

Charge  transfer  coefficients,  aA  and  ac 

1.0  [10] 

Electrolyte  ionic  conductivity  for  ZrP-PTFE, 

5.0  Sm-1 

ctely 

Equilibrium  potential  of  conductor  phase  at 

0.136V  [1] 

the  anode,  <Ppt-EQ,A 

Equilibrium  potential  of  conductor  phase  at 

1.229V 

the  cathode,  ^pt-EQ,c 

Equilibrium  potential  of  electrolyte  phase, 

0.136V 

^ELY-EQ 

Apparent  bulk  density  of  catalyst  support,  pCAT 

0.259  gcatalyst  mLcatalyst 

Specific  surface  area  of  catalyst  support  in  the 

255  m2  <y— 1 

catalyst  "catalyst 

anode  and  cathode,  ACat 

Gas  phase  volume  fraction  in  anode  and 

0.5 

cathode,  eG 

Electrolyte  phase  volume  fraction  in  anode  and 

0.4 

cathode,  £Ely 

Effective  particle  diameter  in  anode  and 

5  pan 

cathode,  Dp 

Land  width,  Lw 

1-6  mm 

Anode  and  cathode  thickness,  ThA,  Thc 

200-400  |Jim 

Membrane  thickness,  ThM 

100  pan 

Fluid  channels  width  in  bi-polar  plates 

0.4  mm 

The  partial  differential  equations  and  the  required  boundary 
conditions  discussed  in  Sections  2.1  and  2.2,  respectively,  are 
rewritten  into  variational  form,  and  discretized  by  the  Finite  Ele¬ 
ment  Method  in  space.  The  Euler  time  stepping  method  is  used  to 
discretize  the  equations  in  time  with  all  nonlinear  reaction  terms 
taken  explicitly  and  linearized  differential  terms  taken  implicitly 
to  increase  stability.  The  partial  differential  equations  are  solved 
one  after  another  in  a  decoupled  fashion.  Time  steps  are  performed 
until  a  steady-state  solution  of  the  fully  coupled  system  of  nonlin¬ 
ear  equations  is  reached,  usually  of  the  order  of  a  few  hundred 
time-steps.  Grid  independence  of  the  solution  was  achieved  by 
computing  solutions  on  sufficiently  refined  meshes. 

A  special  method  was  required  to  couple  the  solutions  on  all  the 
layers.  One  of  the  main  difficulties  is  the  solution  of  the  poten¬ 
tial  equations  for  protons  over  the  three  layers  while  matching 
potentials  and  fluxes  of  protons  at  the  membrane-catalyst  layer 
interfaces.  Moreover,  proton  fluxes  are  zero  on  the  exterior  bound¬ 
ary  of  the  whole  computational  domain  suggesting  that  the  rate  of 
proton  production  in  the  anode  is  equal  to  the  rate  of  proton  con¬ 
sumption  in  the  cathode.  In  an  attempt  to  couple  the  anode  with 
the  membrane,  and  the  membrane  with  the  cathode,  the  anode 
and  the  cathode  are  coupled  together  as  well.  A  domain  decom¬ 
position  method  with  Robin  boundary  conditions  on  the  potential 
has  been  employed  successfully  to  solve  the  potential  equations 
in  the  anode,  membrane  and  cathode.  In  the  Robin  method,  linear 
combinations  of  the  unknown  potentials  and  their  fluxes  are  used 
as  boundary  conditions.  For  example,  Eqs.  (29)  and  (30)  are  used 
as  boundary  conditions  at  the  anode/membrane  interface  to  solve 
potential  equations  in  the  anode  and  membrane  domains: 


£ELY,A°rELY,A 


f  9<^ely,a\ 

V  ^  J 


+  ^ELY,A  =  OML  (  dy  ) 


+  ai  ^ML 


(29) 


f  d@ML 

v  9y 


+  C^2^ML  =  -£ELY,AOTLY,A 


(  9^ely,a\ 

V  *  J 


+  «2^ELY,A 

(30) 


where  ol\  and  a 2  are  relaxation  parameters  related  to  the  anode 
and  membrane  volume,  respectively. 
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Fig.  4.  Polarization  curves  of  direct  propane/oxygen  fuel  cell  using  Pt  anode  and 
cathode:  (a)  experimental  results  [21]  using  Nation  117  at  95  °C,  (b)  experimental 
results  [9]  using  95%  H3P04  at  200 °C,  (c)  experimental  results  [21]  using  H2SO4 
doped  PBI  at  95  °C,  and  (d)  the  present  model  results  for  a  solid  ZrP-PTFE  electrolyte 
at  150°C;  the  base  case-study. 

The  layers  were  solved  sequentially,  with  a  single  iteration  of  the 
coupling  method  at  each  time  step.  A  Newton  method  was  used  to 
force  equal  proton  production/consumption  rates  at  the  anode  and 
cathode. 

2.5.  Model  validation 

The  model  predicts  the  performance  of  a  DPFC  having  inter- 
digitated  flow  fields  and  a  solid  ZrP-PTFE  electrolyte  operating  at 
temperatures  of  at  least  150  °C.  In  order  to  validate  the  model,  two 
sets  of  published  experimental  data  have  been  chosen  for  com¬ 
parison  purposes.  Savadogo  and  Rodriguez  Varela  [21]  reported 
polarization  curves  for  DPFCs  with  several  modified  solid  elec¬ 
trolytes  operating  at  low  temperatures.  Grubb  and  Michalske  [9] 
presented  polarization  curve  data  for  DPFCs  using  phosphoric  acid 
as  the  electrolyte  operating  at  temperatures  as  high  as  200  °C. 

Fig.  4  compares  the  modeling  results  with  the  experimental 
data.  Graph  (a)  in  Fig.  4  shows  the  performance  of  a  DPFC  with 
a  non-modified  Nation  117  solid  electrolyte.  Although  Nation  has 
a  proton  conductivity  approaching  lOSnrr1  [25],  it  has  to  be  used 
at  temperatures  of  80  °C  or  less.  At  elevated  temperatures  both  the 
liquid  phase  water  content  of  Nation  membranes  and  their  perfor¬ 
mance  decline  substantially.  Graph  (b)  in  Fig.  4  shows  Grubb’s  best 
results  obtained  using  a  DPFC  having  a  liquid  phosphoric  acid  elec¬ 
trolyte.  Because  there  are  no  experimental  data  for  direct  propane 
fuel  cells  having  zirconium  phosphate  electrolytes,  a  comparison 
was  made  between  results  from  our  model  having  a  zirconium 
phosphate  electrolyte  and  experimental  data  obtained  with  other 
electrolytes.  The  comparison  in  Fig.  4  indicates  that  the  zirconium 
phosphate  electrolyte  would  have  results  that  are  somewhat  com¬ 
parable  to  the  other  electrolytes.  Although  the  variation  in  results 
among  the  electrolytes  shown  in  Fig.  4  for  direct  propane  fuel  cells 
are  significant,  they  are  minor  when  they  are  compared  to  results 
from  a  hydrogen  fuel  cell. 

The  different  results  obtained  with  the  different  electrolytes 
may  be  related  to  the  various  electrolyte  conductivities.  The  proton 
conductivity  of  phosphoric  acid  electrolyte  depends  on  its  con¬ 
centration  and  on  the  operating  temperature.  Dobos  has  reported 
the  proton  conductivity  of  phosphoric  acid  at  298  K  and  its  tem¬ 
perature  variations  [22].  Proton  conductivity  of  a  95%  H3PO4  at 
200  °C  is  equal  to  35  Sm-1.  The  best  experimental  conductivity 
achieved  for  ZrP  in  our  laboratory  was  approximately  5  S  nrr 1 .  The 
greater  proton  conductivity  values  for  Nation  117  and  95%  FI3PO4 


Fig.  5.  Polarization  curves  of  direct  propane/oxygen  fuel  cell  showing  model  results 
for  catalyst  thickness  of  300  |xm  at  different  land  widths. 

result  in  higher  cell  potentials  at  low  current  densities  compared 
to  the  ZrP-PTFE  electrolyte  modeled  in  this  study.  Flowever,  as 
current  densities  increases,  the  experimental  cell  potential  results 
decrease  more  rapidly  for  Nation  117  and  95%  FI3PO4  compared  to 
the  ZrP-PTFE  model  results.  One  explanation  is  that  both  Nation  and 
95%  H3PO4  have  at  least  a  partial  liquid  film  around  the  Pt/C  cata¬ 
lyst  sites.  This  leads  to  larger  Ohmic  losses  compared  to  a  DPFC  with 
solid  ZrP-PTFE  electrolyte  at  150  °C  where  there  is  no  liquid  phase 
water  and  no  liquid  film  is  present  to  create  a  diffusion  resistance. 

The  model  prediction  in  Fig.  4  indicates  that  a  ZrP-PTFE  mem¬ 
brane  will  give  a  polarization  curve  that  is  generally  similar  to  those 
obtained  with  FI3PO4,  Nation,  and  H2S04  doped  polybenzimidazole 
(PBI).  The  difference  is  that  ZrP-PTFE  does  not  need  liquid  phase 
water,  whereas  the  others  do.  Acceptable  operation  in  the  absence 
of  liquid  water  suggests  that  the  highly  corrosive  environment  can 
be  avoided  and  that  a  search  for  a  non-precious  metal  catalyst  to 
replace  Pt/C  may  be  warranted. 

3.  Results  and  discussion 

The  model  was  employed  for  two  purposes.  The  effects  of  some 
of  the  physical  characteristics  of  the  flowfields  (land  width)  and  the 
MEA  (catalyst  layer  thickness)  were  examined.  It  was  also  used  in  an 
attempt  to  understand  proton  transport  in  the  ZrP  electrolyte  phase 
and  electron  transport  in  the  Pt/C  phase.  In  addition,  the  variation 
of  some  of  the  gas  phase  variables  was  investigated.  These  studies 
are  described  in  this  section. 

Fig.  5  shows  the  effect  of  land  width,  i.e.  the  distance  between 
inlet  and  outlet  channels,  on  the  cell  performance.  It  is  apparent  that 
increasing  the  land  width  has  a  small  negative  effect  on  the  polar¬ 
ization  curve.  As  the  reactants  flow  from  an  inlet  channel  towards 
its  adjacent  outlet  channel  within  the  catalyst  layer,  the  partial 
pressure  of  the  reactants  decreases  and  that  leads  to  a  lower  reac¬ 
tion  rate  and  a  smaller  current  production.  This  can  be  explained  by 
the  dependence  of  exchange  current  density  on  the  partial  pressure 
of  the  reactants  in  the  Butler-Volmer  equation. 

A  major  concern  about  the  interdigitated  flow  fields  is  their 
higher  pressure-drop  compared  to  the  conventional  serpentine 
flow  fields.  Flowever,  it  has  been  shown  that  with  interdigitated 
flow  fields,  thicker  catalyst  layers  have  smaller  pressure  drops  [7]. 
Large  pressure  drops  are  a  serious  drawback  for  interdigitated  flow 
fields.  This  perception  is  correct  for  serpentine  flow  fields  using 
conventional  anode  catalyst  layers  of  40  pum.  Flowever,  by  increas¬ 
ing  the  thickness  of  the  anode  catalyst  layer  to  400  p,m,  the  pressure 
drop  along  the  land  width  will  decrease.  Moreover,  an  increase  in 
the  catalyst  thickness  results  in  an  improvement  in  the  polariza- 
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Fig.  6.  Polarization  curves  of  direct  propane/oxygen  fuel  cell  showing  model  results 
for  the  land  width  of  4  mm  at  different  catalyst  thicknesses. 


Fig.  7.  Potential  profile  for  the  ZrP  phase  within  the  cathode  catalyst  layer,  the  anode 
catalyst  layer,  and  the  membrane  layer. 


tion  curve  as  can  be  seen  in  Fig.  6.  Improved  dispersion  of  the  metal 
catalyst  should  be  obtained  using  a  thicker  catalyst  layer  without 
changing  the  total  amount  of  catalyst  used.  An  increased  dispersion 
would  correspond  to  a  larger  number  of  reaction  sites. 

Fig.  7  shows  the  two-dimensional  variation  of  the  electrical 
potential  in  the  ZrP  phase  of  the  entire  domain,  that  is  the  anode 
catalyst  layer,  the  membrane  layer,  and  the  cathode  catalyst  layer. 
A  cross-section  of  Fig.  7  is  shown  in  Fig.  8,  where  the  electrical 
potential  in  the  ZrP  phase  is  compared  with  that  in  the  Pt/C  phase. 

Fig.  8  shows  the  following  MEA  profiles  in  the  y-direction  (as 
defined  in  Fig.  2):  (a)  a  comparison  of  the  electrical  potential  in  the 
Pt/C  solid  phase  with  that  in  the  ZrP  phase  for  the  anode  catalyst 
layer,  (b)  the  same  comparison  as  in  (a)  but  for  the  cathode  catalyst 
layer,  and  (c)  the  electrical  potential  profile  of  the  ZrP  phase  in  the 
membrane  layer.  As  the  electrical  conductivity  in  the  Pt/C  phases 
is  very  high,  the  electrical  potential  in  both  of  the  Pt/C  phases  in 


- Electrolyte  phase  - Conductor  phase 


y-axis  /  mm 


Fig.  8.  Potential  profiles  in  the  y-direction  for  the  ZrP  and  Pt/C  phases  located  at 
the  middle  of  the  domain  x-direction  for  the  cathode  and  anode  catalyst  layers  and 
membrane  layer.  The  arrows  point  in  the  direction  of  the  ordinate  scale  that  applies 
to  each  of  the  three  curves. 


their  respective  anode  and  cathode  layers  are  almost  constant,  as 
shown  by  the  horizontal  dashed  lines  in  Fig.  8.  However,  the  electri¬ 
cal  potential  in  the  ZrP  phase  varies  according  to  Eqs.  (6)  and  (14), 
as  shown  by  the  solid  line  in  Fig.  8. 

The  model  uses  this  electrical  potential  gradient  to  account 
for  proton  migration  from  the  ZrP  phase  in  the  anode  catalyst 
layer,  through  the  ZrP  phase  in  the  membrane  layer,  and  into  the 
ZrP  phase  in  the  cathode  catalyst  layer.  As  mentioned  previously 
neglecting  the  diffusion  term  in  Eqs.  (6)  and  (14)  is  a  common 
assumption  used  in  most  PEM  fuel  cell  models.  Reasonable  pre¬ 
dictions  of  the  fuel  cell  polarization  curve  were  obtained  using  this 
assumption,  as  shown  in  Figs.  4-6.  However,  one  of  the  undesir¬ 
able  consequences  of  this  assumption  is  that  the  slope  of  electrical 
potential  in  the  ZrP  phases  in  the  three  layers  is  positive.  It  should 
be  negative.  The  correct  sign  of  the  slope  can  be  understood  from 
the  following  reasoning:  the  potential  gradient  between  anode  and 
cathode  drives  electrons  from  the  Pt/C  phase  of  the  anode  to  the 
Pt/C  phase  of  the  cathode  via  the  external  circuit.  If  the  fuel  cell 
had  an  aqueous  electrolyte,  the  electrical  double  layer  would  have 
described  the  difference  in  electrical  potential  between  the  catalyst 
layers  and  the  electrolyte.  Since  ZrP  has  semi-conducting  properties 
[23],  the  difference  in  electrical  potential  between  Pt/C  and  ZrP  will 
be  described  by  a  space  charge  region  [24].  As  a  result  the  electrical 
potential  in  the  ZrP  phases  will  be  directly  related  to  the  electrical 
potentials  in  the  Pt/C  phases  that  are  shown  as  the  dashed  lines  in 
Fig.  8.  Therefore  the  electrical  potential  in  the  ZrP  phase  of  the  anode 
catalyst  layer  must  be  less  than  that  in  the  ZrP  phase  of  the  cath¬ 
ode  catalyst  layer.  Therefore  opposite  trend,  shown  in  Fig.  8,  cannot 
be  correct.  Since  negatively  charged  electrons  are  driven  from  the 
anode  to  the  cathode  by  the  electrical  potential  gradient  (in  the 
Pt/C  phases),  the  corresponding  gradient  (in  the  ZrP  phases),  when 
correctly  calculated,  must  provide  a  driving  force  for  positively 
charged  protons  from  the  cathode  to  the  anode.  As  the  electrochem¬ 
ical  reaction  proceeds  protons  will  accumulate  (in  the  ZrP  phase)  at 
the  anode  until  the  proton  concentration  gradient  from  the  anode 
to  the  cathode  is  sufficient  to  overcome  the  electrical  potential  gra¬ 
dient  (in  the  ZrP  phases).  When  that  condition  is  attained,  electrons 
will  flow  through  the  ZrP  phases  from  the  anode  to  the  cathode.  The 
conclusion  is  that  Fig.  8  has  an  incorrect  sign  for  the  slope  of  the 
electrical  potential  gradient  (in  the  ZrP  phase).  If  the  model  had 
included  the  diffusion  term  with  its  concentration  gradient  in  Eq. 
(21 )  then  the  sign  of  the  slope  in  Fig.  8  would  have  been  correct. 

It  is  worth  mentioning  that  the  solid  line  in  Fig.  8  for  the  elec¬ 
trolyte  phase  potential  shows  that  the  condition  of  zero  proton 
flux  at  the  lands  has  been  satisfied,  and  coupling  at  the  mem¬ 
brane/catalyst  layers  interfaces  has  been  done  perfectly. 

Fig.  9  shows  that  the  overpotential  predicted  by  the  model 
at  both  the  membrane/catalyst  interfaces  is  greater  than  that 
at  their  corresponding  lands.  This  prediction  by  the  model  is 
incorrect,  as  can  be  seen  from  the  following  reasoning  at  the 
anode:  the  overpotential  at  the  anode  is  given  by  Eq.  (9), 
rjA  =  (^pt,A  -  ^ely,a)  -  (^Pt-EQ,A  -  ^ely-eq.)  =  (^pt,A  -  ^ely.a)  when 
the  equilibrium  potential  of  the  ZrP  electrolyte,  ^ely-eq  is  been 
defined  to  be  equal  to  the  equilibrium  potential  in  the  Pt/C  catalyst 
^pt-EQA-  According  to  Fig.  8,  the  potential  difference  0ptA  -  ^ely,a 
is  greater  at  the  anode  catalyst-membrane  interface  than  at  the 
anode  catalyst-land  interface.  That  is  consistent  with  the  overpo¬ 
tential  in  Fig.  9.  However,  if  proton  diffusion  driven  by  the  proton 
concentration  gradient  had  been  included  in  the  model,  then  the 
electrical  potential  gradient  in  the  ZrP  phases  in  Fig.  8  would  have 
been  reversed.  That  would  also  reverse  the  slope  of  the  overpoten¬ 
tial  gradients  in  Fig.  9. 

Fig.  10a  and  b  shows  the  magnitude  and  direction  of  proton  and 
electron  fluxes,  respectively.  Protons  and  electrons  are  produced 
in  the  anode  and  consumed  in  the  cathode.  Protons  travel  through 
the  ZrP-PTFE  electrolyte  phase  from  anode  to  cathode,  as  shown  in 
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y-axis  /  mm 


Fig.  9.  Overpotential  profile  at  the  anode  and  cathode. 


Oxygen  mole  fraction 
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Fig.  11.  Propane  and  oxygen  mole  fraction  at  the  anode  and  cathode,  respectively. 


Fig.  10a.  Protons  are  produced  throughout  the  anode  catalyst  layer. 
The  arrow’s  length  in  Fig.  10a  correctly  indicates  that  the  proton 
flux  increases  as  the  position  in  the  anode’s  y-direction  changes 
from  the  anode  catalyst-land  interface  to  the  membrane-anode 
catalyst  interface. 

Electrons  travel  from  anode  to  cathode  via  the  external  circuit 
although  there  is  a  very  small  amount  of  electronic  current  through 
the  membrane  layer.  Electron  transport  also  occurs  in  the  ZrP  phase 
of  the  catalyst  layers.  It  is  the  electron  transport  in  the  ZrP  phase 
that  is  shown  in  Fig.  1 0b.  Neglecting  the  diffusion  term  in  the  proton 
conservation  equation  caused  the  direction  of  the  electron  flux,  in 
Fig.  10b,  to  be  incorrect,  as  was  described  previously. 

Fig.  1 1  shows  the  variation  in  the  propane  and  oxygen  mole  frac¬ 
tions  at  the  anode  and  cathode,  respectively.  The  concentrations  of 
the  reactants  decrease  along  the  land  as  the  reactants  flow  from 
the  inlet  to  the  outlet.  Therefore,  changing  the  land’s  width  can 
control  the  conversion  of  the  reactants.  This  is  an  important  bene¬ 
fit  of  interdigitated  flow  fields  for  DPFCs.  Even  if  the  feed  channels 
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Fig.  10.  (a)  Protonic  current  from  anode  to  cathode  in  electrolyte  phase.  The  vectors 
length  indicates  the  current  magnitude  which  varies  from  0  to  120  mAcm-2  in  this 
case  and  (b)  electronic  current  from  cathode  to  anode  in  electrolyte  phase.  The  vec¬ 
tors  length  indicates  the  current  magnitude  which  varies  from  0  to  2e  -  14  mAcm-2 
in  the  same  case  as  Fig.  10a. 


Fig.  12.  Velocity  vectors  showing  movement  of  fluid  in  the  anode  and  cathode 
catalyst  layers. 

are  dead-ended  so  that  all  of  the  reactants  must  flow  through  the 
catalyst  layer,  as  shown  in  Fig.  1,  it  would  be  undesirable  for  unre¬ 
acted  propane  to  contaminate  the  product  gas  stream  that  ideally 
should  consist  only  of  carbon  dioxide  and  water  vapour.  The  water 
concentration  at  the  anode  entrance  was  specified  to  exceed  the 
stoichiometric  ratio  because  it  is  needed  for  both  reaction  and  to 
assist  proton  transport  (osmotic  drag). 

Fig.  12  shows  the  velocity  vectors  in  the  gas  phase  of  the  anode 
and  cathode  catalyst  layers.  Arrows  in  x-direction  on  the  land  rep¬ 
resent  a  slip  flow  regime  in  the  catalyst  layer.  This  suggests  that 
the  Ergun  equation  describes  momentum  conservation  appropri¬ 
ately,  since  the  catalyst  layers  have  been  modeled  as  packed  beds 


Fig.  13.  Propane  conversion  at  different  anode  and  cathode  land  widths. 
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and  flow  though  a  packed  bed  is  possible  even  at  the  wall  of  a  land. 
Finally,  the  small  velocity  vectors  indicate  poor  catalyst  turnover 
frequency. 

Fig.  13  shows  that  the  propane  conversion  increases  as  the 
land  width  increases.  The  fact  that  100%  propane  conversion  can 
be  approached  with  interdigitated  flow  fields  is  important.  As 
explained  previously  unreacted  propane  in  the  product  stream 
would  be  undesirable.  A  100%  conversion  result  is  unlikely  with  a 
serpentine  flow  field.  In  a  serpentine  flow  field  the  carbon  dioxide 
products  would  mix  with  the  hydrocarbon  feedstock.  The  smaller 
propane  concentration  in  the  channels  of  serpentine  flow  fields 
would  decrease  the  driving  force  that  pushes  the  hydrocarbons 
into  the  catalyst  layer. 

4.  Conclusion 

The  performance  predicted  for  a  DPFC  at  150°C  using  a  polyte- 
trafluoroethylene  membrane  with  its  pores  filled  with  zirconium 
phosphate  (ZrP-PTFE)  was  similar  to  experimental  results  obtained 
with  Nation,  aqueous  FI3PO4,  and  FI2SO4  doped  polybenzimidazole 
electrolytes,  when  they  were  used  in  DPFCs.  Because  a  ZrP-PTFE 
membrane  operating  at  1 50  °C  does  not  require  water  in  the  liquid 
phase,  corrosion  will  be  much  less  severe,  and  the  use  of  non¬ 
precious  metal  catalysts  may  be  possible.  It  was  found  that  the 
thickness  of  the  catalyst  layer  could  be  increased  sufficiently  so 
that  the  pressure  drop  between  the  reactant  and  product  channels 
of  interdigitated  flow  fields  was  small.  By  increasing  the  width  of 
the  lands  and  therefore  the  reactant’s  contact  time  with  the  cata¬ 
lyst  it  was  possible  to  approach  100%  propane  conversion,  thereby 
ensuring  that  there  will  be  a  minimum  concentration  of  propane 
in  the  product  stream  from  the  fuel  cell.  Finally  neglecting  the  dif¬ 
fusion  term  in  the  equation  for  the  conservation  of  protons  causes 
serious  errors  in  computations  of  the  electrical  potential  in  the  ZrP 
phase,  the  electron  flux  in  the  Pt/C  phase,  and  the  overpotential  in 
both  the  anode  and  cathode  catalyst  layers. 
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